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Abstract 

The coincidence problem of late cosmic acceleration is a serious riddle in connection with our 
understanding of the evolution of the Universe. In this paper we show that an interaction between 
the dark energy component (either phantom or quintessence) and dark matter can alleviate it. In 
this scenario the baryon component is independently conserved. This generalizes a previous study 
[S. del Campo, R. Herrera, and D. Pavon, Phys. Rev. D 71, 123529 (2005)] in which neither 
baryons nor phantom energy were considered. 
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I. INTRODUCTION 



According to the conventional picture, the present accelerated expansion of the Universe 
is driven by the negative pressure of an unknown and unclustered component (dubbed 
"dark energy") that currently contributes about 70% of the total density. The remaining 
30 per cent is shared between cold dark matter (pd m ~ 25%) and cold baryons (pb ~ 
5%) Q fl . The two latter, being pressureless, redshift with expansion faster than 

the dark energy component. Then, the cosmic coincidence problem arises, "Why are the 
densities of matter and dark energy of precisely the same order today?" p. Clearly, this 
is one outstanding riddle in our understanding of the Universe. To solve it one is forced to 
adopt an evolving dark energy field (either quintessence or phantom energy) or accept an 
incredibly tiny cosmological constant and admit that the "coincidence" is just a coincidence 
that hopefully might be somewhat alleviated with the help of the anthropic idea Q|. Here 
we take the view that before resorting to the second option we should further explore the 
first one. Yet, an evolving dark energy cannot solve the problem either unless a suitable 
interaction (coupling) with matter is allowed 0,0]. Note that the coupling alters the rate at 
which both matter and dark energy redshift with expansion; this is why it can potentially 
alleviate the aforesaid problem. 

Interestingly, rather than in connection with the coincidence problem, which was not even 
formulated at the moment, this interaction was first proposed as a mechanism to reduce the 
value of the cosmological constant [loj ]. Notice that in the absence of underlying symmetry 
that would suppress the coupling matter-dark energy there is no a priori reason to dismiss 
it. In the last years, various proposals at the fundamental level for the coupling leading 
to a constant ratio matter/dark energy at late times were advanced (see e.g., ll| and 
references therein), and specific phenomenological models have been built and contrasted 
with observational data (high redshift supernovae and CMB anisotropies) and seen to pass 
the tests [l^. Further, the Akaike and Bayesian informative criteria ^| when applied 
to high redshift supernovae data suggest a transfer of energy from the phantom component 
to the matter component; yet, the conventional ACDM model is still preferred 

While a constant ratio matter /dark-energy at late times (including the present one) would 
clearly alleviate the coincidence problem it should be noted that a much less strong condition 
would serve. It would suffice that nowadays the aforesaid ratio varies only slowly (i.e., less 
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faster than the scale factor) and be of order unity ( "soft coincidence" ) . Recently, it was found 
within this approach that the quintessence scenario was favored over the tachyon scenario 
for the latter would imply an excessive amount of pressureless matter today [8|. However, 
in order to circumvent the tight constraints on long-range forces [16] the baryon component 
was left aside altogether. This may be justified because the analysis was restricted to times 
about the present one which, as said above, is characterized, among other things, by a low 
value to the baryon density. In this paper we generalize the analysis by including baryons in 
the energy density budget as an independently conserved component and extend the study 
to earlier times -though not to the radiation dominated era. As dark energy components 
quintessence and phantom are separately considered. 

The outline of the paper is as follows. In section II we present the model. In section III we 
constrain it with recent high redshift supernovae data. Finally, in section IV we discuss and 
summarize our findings. As usual, a zero subscript or superscript attached to any quantity 
means that it should be evaluated at the present epoch. 

II. THE INTERACTING MODEL 

We consider a spatially flat Friedmann-Lemaitre-Robertson- Walker universe dominated 
by a three-component system, namely, cold baryonic matter, cold non-baryonic dark matter 
and dark energy, such that the two latter components do not conserve separately but interact 
with each other in a manner to be specified below. The energy density and pressure of the 
dark energy, assuming it is a quintessence field, are given by 

P<t> = \^ + V(<P), and P^ = l0 2 -K(0), (1) 

respectively. If the dark energy is a phantom field we have instead, 

P4> = -±0 2 + V(<j>) , and p^-l^-V^). (2) 

The upper-dot stands for derivative with respect to the cosmic time and V(4>) denotes both 
the quintessence field potential and phantom potential. As is usually done, we postulate that 
the dark energy component (either quintessence or phantom) obeys a barotropic equation of 
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state, i.e., P^ = w^p^ with a negative constant of order unity (a distinguishing feature 
of dark energy fields is a high negative pressure). 

We assume that the dark matter and dark energy components are coupled through a 
source (loss) term (say, Q) that enters the energy balances 

Pdm + 3Hp dm = Q, (3) 



and 



+ 3H{p^ + P^ 



-Q 



(4) 



In view of Eq.Q last expression can alternatively be written as 



+ 3H(f) + V 



-Q. 



In the case that the scalar field is of phantom type the corresponding equation reads 



+ 3H(f) - V 



Q, where the prime denotes derivative with respect to 



We consider that the baryon component is conserved whence its energy density redshifts 



as 



Pb = P b 



ofaoV 



(5) 



Defining p m = Pb + Pdm and using Eqs.(jSJ) and (JHJ) we obtain 



p m + 3Hp m = Q. 



(6) 



Likewise, we assume that the interaction is related to the total energy density of matter 
(baryonic plus dark) by Q = 3c 2 Hp m , where c 2 is a small positive-definite constant. As we 
shall see this choice ensures that the ratio between the energy densities, r(a) = Pm/p^, is 
a monotonous decreasing function of the scale factor, and such that around present time 
it varies very slowly. By very slowly we mean that | (r/r) |^ H Q . This contrasts with 
previous studies in which r was demanded to asymptotically approach a fixed value at late 
times. 

Clearly, 



r 



d 
It 



Pr, 



Pr, 



(7) 
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In virtue of the above expressions last equation boils down to 



r = 3Hr [c 2 (l + r) — , 



(8) 



whose solution reads 




-i -i 



r = r £ c 2 r — 



(9) 



with £ = — c 2 > 0. Figure ^ depicts the monotonous decrease of r from values higher 



than unity at early times to a nearly constant value at present time (which we have fixed 
as 3/7) irrespective of whether the dark energy is a quintessence field or a phantom field. 
Notice that for vanishing scale factor r tends to the finite, constant value £/c 2 . However, 
our model should not be extrapolated to such early stage. 

One may wonder about the size of c 2 . Obviously it should not be large. In any case 
it must be lower than |u^|/(l + r ); otherwise, by Eq. (|SJ), the Universe would have been 
dominated by the dark energy from the beginning of the expansion and galaxies would not 
have come into existence. On the other hand, it should not be very small for it would have 
a negligible impact and our model would be hardly distinguishable (depending on the value 
of w^) from the standard quintessence or phantom models. 

Because of c 2 ought to be a small quantity -though not very small-, the right hand side 
of Eq.© always stays above zero, becoming negligible only for a 3> a . The closer c 2 
is to | |/(1 + r o), the more alleviated the coincidence problem gets. In particular, for 
c 2 > [| We/, | — 1/3] (1 + ro)^ 1 the current rate, | (r/r)o |, results lower than Hq (bear in mind 
that the corresponding rate in case the dark energy were just the cosmological constant is 
3H ), whereby the criterion of "soft coincidence" is satisfied and the coincidence problem 
gets significantly alleviated. 

From Eq.© along with the expression for Q we get 




(10) 
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FIG. 1: Evolution of the ratio Pm/Pcp as given by Eq.® for c 2 = 10~ 5 (upper panel) and c 2 = 10~ 3 
(lower panel). In both panels tq = 3/7 with = —0.8 (quintessence) and w$ = —1.1 (phantom). 
Not shown is the behavior of r for a — > as the model does not apply to such early times. 

Thus, the densities of dark matter and dark energy are given by 




and 



c 2 r 



35 

-) [ c 2 (l + r )-K|] 



(12) 



respectively. Obviously, the constant c 2 could be determined if either p m , given by Eq. (fTUj) 
above, or the ratio Pb/ p m oc a~ 3c2 , were accurately known at different redshifts. The above 
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expression for pd m may suggest that this quantity becomes negative at small scale factor. 
That is so; however, if one takes into account that p\j p m * s about 0.2, for reasonable values 
of c 2 , this only occur well in the radiation era, i.e., beyond the range of applicability of our 
model. 

As Fig. n shows, r > 10 at sufficiently early times (e.g., at redshifts larger than, say, 3). 
This is consistent with the analysis of Caldwell et al. |l7( who found that at the epoch of 
last scattering (z — 1, 100) as well as at the onset of structure formation (z ~ 10 3 ) (the 
energy density of the dark energy in units of the critical density) should not exceed 0.1. 

On the other hand, from Eqs. (JHJ) and © alongside the condition Q > it follows that 
the baryon density never dominates the matter density (i.e., Pb/ Pdm < 1 always). This is in 
keeping with the widely accepted scenario of cosmic structure formation, in which after the 
last scattering the baryonic matter freely fells into the deep potential wells created by the 
dark matter. If p^, were larger than pd m at early times, then the aforesaid scenario could be 
spoiled. 

By combining Friedmann's equation 



SH 2 = K (p b 



(k = 8%G), 



(13) 



with Eqs. (fTU jl - lfT^ we get the Hubble function 



H{a) 



Ho 



V(l+r )f V a 



3? 

KK-(-) [c 2 (l+r )-K|] 

do 



1/2 



(14) 



where H = y (1 + tq)k p°/3, which will be needed below both to obtain the luminosity 
distance -a previous step to draw the likelihood contours- and the deceleration parameter, 
q = —a/(aH 2 ). The evolution of the latter as a function of the redshift is depicted in Fig. 

121 

In this model the universe began accelerating only recently, at redshifts about 0.65. As we 
have checked numerically this behavior is scarcely sensitive to the value of the c 2 parameter 
provided this lies in the range 10 -5 < c 2 < 10 _1 (this is why we content ourselves with 
presenting just two plots of q vs. z). The larger c 2 , the larger the transition redshift. 

Using Eq.([IJ) (respectively, Eq.(j2J) the potential for the quintessence field (respectively, 
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FIG. 2: The deceleration parameter as a function of the redshift, z = {clq/o) — 1. Here ro = 3/7 
and c 2 = 1CT 3 . 

phantom field) in terms of the scale factor reads 



V{<j>{a)) = (1 + K|) 



a 



(15) 



2r ^ V a 

To obtain the potential as a function of we must first express the latter in terms of a. 
To this end, we write 



L2_ ± , z^Ha 
da 



(16) 



where the plus(minus) sign corresponds to quintessence field (phantom field). And in virtue 
of Eqs. ([T2 |) . (|T3 |) and (fT51) it follows that 



±3(1 -K|) 1 



3£ 



In 3(a) + 



c 



(17) 



S 



where 

3(a) = -cV + 2c 2 (a/a ) 35 (l + r ) - (r + 2(a/a ) 3f )K| 
+2 A /|^|(a/a ) 3 « - C 2((a/a )3« + r [(a/a ) 3 € - 1]) ^/|^|(r + (a/a ) 3 ?) - c 2 (l + r )(a/a ) 3 ? 



a) = 3 f a/ !^^ [ c 4 (l + r )(a/a ) 3 « + c 2 r K| [(a/a ) 3 « - 2] - ^(a/a ) 



3£ 



C 3 



2c^\^\w^\{a/a^ - c 2 [(a/a ) 3 « + r ((a/a ) 3 « - 1)] 



x ^ (r + (a/a o ) 3? )|w | - c 2 (l + r )(a/a ) 3 ^ 



and C denotes the integration constant 



c 2 



C = \n%(a ) + J- lnK(ao). 

Kl 

The dependence of the potential on the field <fi is shown in Fig. El 

Clearly, this potential must be understood as an effective one. It should be noted that 
these effective potentials (quintessence and phantom potentials) are similar to those used in 
inflationary models, where the universe undergoes an accelerated period at very early times. 
Thus, our potentials describe an accelerated phase at present time, and simultaneously 
alleviate the coincidence problem. 

III. COMPARING WITH SUPERNOVAE DATA 

In this section we use two independent supernovae type la (SNIa) data sets, namely 
the "gold" sample compiled by the High-Z Supernovae Search Team (HZT) [3], and the 
Supernova Legacy Survey (SNLS) sample [4j, to constrain the parameters of the model. 
The "gold" sample, collected from different sources, comprises 157 SNIa, of redshifts up to 
1.5 (14 of which, discovered by the Hubble Space Telescope, lie in the interval 1 < z < 1.5), 
with reduced calibration errors coming from systematics. The SNLS sample is smaller, 71 
SNIa, with redshifts below unity. However, the technique employed guarantees that no 
source is lost and the data are of a higher quality. 
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FIG. 3: The normalized potential, V{4>)/p^ n , vs y/K<p. Here tq = 3/7, c 2 = 10 3 and 0o = 
(C/3£) y/ ±3(1 — {w^D/k, the plus (minus) sign corresponds to the quintessence field (phantom 
field). 

Figures |U and 03 depict the best-fit of the phantom model (solid line) and the quintessence 
model (dashed line) to the "gold" sample and the SNLS sample, respectively. For the sake 
of comparison the flat ACDM model is also shown. In plotting the graphs the expression 
for the distance modulus, \i — 5 log d L + 25, was employed. Here di = (1 + z) f* H~ x {z'}dz' ', 
is the luminosity distance in megaparsecs. 

Figures El and [7| portray the two-dimensional likelihood contours for the case that the 
dark energy component is a quintessence field, based on the "gold" P and the SNLS sample 
jj], respectively. Both set of contours were calculated by running through a grid of models 
on a four- dimensional parameter space. The prior Q m + = 1 was used and the present 
value of the Hubble parameter was allow to vary in the interval 60 < Hq < 70km/s/Mpc. 
The rest of the priors are 10~ 5 < c 2 < 0.3, and — 1 < w<f, < —0.6. The constraints obtained 
on the free parameters from the HZT data are: = 0.74+^, w^ila) < —0.9, c 2 (lcr) < 
0.13. The latter parameter shows large degeneracy and only in this case we find an upper 
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FIG. 4: Distance modulus vs redshift for the best fit quintessence model, = 0.71, = —0.99, 
c 2 = 10~ 5 , x 2 = 177 (dashed line), and the phantom model, = 0.40, = -4.63, c 2 = 0.20, 
X 2 = 173 (solid line). The flat ACDM model is also included, ft A = 0.70, w A = -1, x 2 = 178 
(dotted line). The data points correspond to the "gold" sample of SNIa 

limit. In its turn, the constraints obtained on the free parameters from the SNLS data 
are: = 0.83+q;qq, w^la) < —0.62. Unfortunately, the c 2 parameter is practically not 
constrained in this case. 

Similarly, Figs. |H1 and E] show the corresponding contours for the phantom model using 
identical sets of data and the priors 1CT 5 < c 2 < 0.3, — 5 < < —0.6. The constraints 
obtained on the free parameters from the first set (HZT) are: = 0.48 ± 0.06, w$ = 
—3.0 ± 1.1. Again, the interaction parameter c 2 shows a wide degeneracy. In its turn, the 
constraints obtained from the SNLS data set are: = 0.71+qj 2 , = —1.1 ± 0.4, with no 
real constraint on c 2 given its large degeneracy. 

As is usual in scenarios of late acceleration, this model prefers a lower contribution of 
dark energy. The fact that phantom models prefer a lower value of is only natural given 
their low value of w^. It interesting to see that the SNLS data favor a higher value for 
(both for phantom and quintessence models) than HZT's. 
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FIG. 5: Distance modulus vs redshift for the best fit quintessence model, fi^ = 0.74, uu = —0.99, 
c 2 = 10~ 5 , x 2 = 111-09 (dashed line), and the phantom model, = 0.70, = -1.10, c 2 = 0.035, 
X 2 = 111.02 (solid line). The flat ACDM model is also included, Q A = 0.72, w A = -1, X 2 = 111.03 
(dotted line). The data points correspond to the SNLS sample Q|. 

IV. DISCUSSION 

We studied a model of late cosmic acceleration by assuming that the dark matter and 
dark energy components are coupled to each other so that there is a transfer of energy from 
the latter to the former, while the baryon component is conserved. By suitably choosing 
the interaction term, Q, the ratio between both dark sources of gravity is seen to evolve, at 
present time, less faster than the scale factor. This considerably alleviates the coincidence 
problem albeit it does not solve it in full. Clearly, to achieve the latter one should derive 
the present value of the aforesaid ratio, or at least show that it has to be of order unity. 
For the time being, r$ ought to be understood as an input parameter. This also holds for a 
handful of key observational quantities such as the present value of the cosmic background 
radiation temperature, the cosmological constant H Q , or the ratio between the number of 
baryons and photons. 
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FIG. 6: Likelihood contours for the quintessence model displaying the 68%, 95%, and 99.99% 
confidence intervals. The likelihood are marginalized over the rest of the parameters. The bottom 
panels show the the probability functions for the quintessence energy density normalized to the 
critical density (left panel), and the equation of state parameter of the quintessence fluid (right 
panel). The data points correspond to the "gold" sample of SNIa 

The transition deceleration-acceleration occurs recently, at redshifs lower than unity (see 
Fig. 12). This contrasts with other interacting models in which the transition is predicted to 

n 

take place much earlier, at redshifhts as hig as 10 |18| . 

It is noteworthy that the expression for the potential, Eq. (|15J). is identical irrespective 
of whether the dark energy component is a phantom or a quintessence field. 

Both when the dark energy is a quintessence or a phantom field the model fits rather well 
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FIG. 7: Same as Fig. El except that the data points correspond to the SNLS data set of Ref. 



the HZT data set, and in the latter case clearly better than the concordance model ACDM 
does (y£w om = 173, ^ quintessence = 177, x\ = 178). Yet, the Bayesian information criterion 
(BIC) given by the formula BIC= x 2 +P m N, where p is the number of free parameters 
of the model (2 for the ACDM model, 4 for dark energy models) and N the number of data 
points, distinctly favors the ACDM model for it yields a lower figure (BIC gu j n f essence ~ 197, 
BlC phantom ~ 193, BIC A ~ 188). 



The fits to the SNLS data are rather similar (xt 



phantom 



111.02, x 



quintessence 



111.09, 



Xa 



111.03). However, the Bayesian information criterion again sides with the ACDM 



model (BICquint essence ~ BIC p w om ^ 128, BIC A ~ 119). 

We have not used the Akaike criterion because it favors models with larger number of 
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FIG. 8: Likelihood contours for the phantom model displaying the 68%, 95%, and 99.99% con- 
fidence intervals. The likelihood are marginalized over the rest of the parameters. The bottom 
panels show the the probability functions for the phantom energy density normalized to the critical 
density (left panel), and the equation of state parameter of the phantom fluid (right panel). The 
data points correspond to the "gold" sample of SNIa k|. 



parameters 




Thus, from the point of view of the Bayesian information criterion the ACDM model is 
preferred, and the soft coincidence model for quintessence seems clearly disfavored if not 
discarded. As for the phantom model, the situation appears somewhat undecided. While, 
regarding the BIC, it lags behind the concordance model it shows a better fit to the HZT 
sample of supernovae and alleviates the coincidence problem. We may summarize by saying 
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that the SNIa data are not yet conclusive thereby a richer supernovae statistics, especially 
at redshifts larger than unity, is needed to come to a verdict. 

Unfortunately, the supernovae constraints on the c 2 parameter is rather poor. This is 
coherent with the fact that the interaction dark energy-dark matter affects the luminosity 
distance only at third order in the redshift j^J. Nevertheless, we hope to be able to break 
the degeneracy by submitting the model to further tests such as the CMB temperature 
anisotropies and the distribution of matter at cosmological scales. This will be the subject 
of a future work. 
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